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Abstract. In this paper we review some recent results concerning inverse 
C/} . problems for thin elastic plates. The plate is assumed to be made by non- 

homogeneous linearly elastic material belonging to a general class of anisotropy. 
A first group of results concerns uniqueness and stability for the determina- 
tion of unknown boundaries, including the cases of cavities and rigid inclu- 
sions. In the second group of results, we consider upper and lower estimates 
of the area of unknown inclusions given in terms of the work exerted by a 
couple field applied at the boundary of the plate. In particular, we extend 
previous size estimates for elastic inclusions to the case of cavities and rigid 
inclusions. 
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1 Introduction 

The problems we consider in the present paper belong to a more general issue 
that has evolved in the last fifteen years in the field of inverse problems. 
Such an issue collects the problems of determining, by a finite number of 
boundary measurements, unknown boundaries and inclusions entering the 
boundary value problems for partial differential equations and systems of 
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elliptic and parabolic type. Such problems arise in nondestructive techniques 
by electrostatic measurements |Ka-Saj . [In] , in thermal imaging [Vo-Moj . 
[EEC], in elasticity theory [N], [Bo-Co] . [M-R3] . and in many other similar 
applications [Is2] . In this paper we try to enlighten the different facets of the 
issue fixing our attention on the theory of thin elastic plates. In Section [3] we 
give a self contained derivation of the Kirchhoff-Love plate model on which 
such a theory is based. 

We begin with the problem of the determination of a rigid inclusion em- 
bedded in a thin elastic plate. 

Let Q denote the middle plane of the plate. We assume that Q is 
a bounded domain of IR 2 of class C 1 ' 1 . Let h be its constant thickness, 
h « diam(fi). The rigid inclusion D is modeled as an open simply con- 
nected domain compactly contained in Q, with boundary of class C 1 ' 1 . The 
transversal displacement w G H 2 (Q) of the plate satisfies the following mixed 
boundary value problem, see, for example, [Fi] and [Gu] . 

' div(div(PV 2 w)) = 0, in Q \ D, q-q 

(PV 2 w)n ■ n = -M n , on dQ, (1.2) 

< div(PV 2 w) -n+ ((PV 2 w)n-r), s = (M T ), S , on dQ, (1.3) 

w\pEA, inD, (1.4) 

w % - w ]m 011 dD, (1.5) 

coupled with the equilibrium conditions for the rigid inclusion D 
I (div(PV 2 w) • n + ((PV 2 M))n • t), s ) g - ((PV 2 w)fi • n)g. n = 0, 

JdD 

for every g G A, (1.6) 

where A denotes the space of affine functions. In the above equations, n and 
r are the unit outer normal and the unit tangent vector to Q\D, respectively, 
and we have defined w e = wlnyo and w l = w\-jy. Moreover, M T , M n are the 

twisting and bending components of the assigned couple field M, respectively. 
The plate tensor P is given by P = f^C, where C is the elasticity tensor 
describing the response of the material of the plate. We assume that C 
has cartesian components Cy^, i,j,k,l = 1,2, which satisfy the standard 
symmetry conditions (14.21) . the regularity assumption (14.31) and the strong 
convexity condition (14.41) . 

Given any M G H~^(dQ,M. 2 ), satisfying the compatibility conditions 
J m Mi = 0, for i — 1,2, problem (ll.ll) - (ll.6p admits a solution w G H 2 (Q), 
which is uniquely determined up to addition of an affine function. 
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Let us denote by F an open portion within dfl representing the part of 
the boundary where measurements are taken. 

The inverse problem consists in determining D from the measurement 
of w and w tU on T. For instance, the uniqueness issue can be formulated as 
follows: Given two solutions Wi to (ll.lB - fll.6p for D = D i} i — 1, 2, satisfying 

wi = w 2 , on T, (1.7) 

wi, n = w 2 , n , on T, (1.8) 

does D\ = D 2 hold? 

It is convenient to replace each solution Wi introduced above with Vi = 
Wi — gi, where <?j is the affine function which coincides with Wi on dDi, 
i = 1,2. By this approach, maintaining the same letter to denote the solution, 
we rephrase the equilibrium problem fll.ip — ()1.5p in terms of the following 
mixed boundary value problem with homogeneous Dirichlet conditions on 
the boundary of the rigid inclusion 

div(div(PV 2 w)) = 0, in Q \ D, (1.9) 

(PV 2 w)n • n = -M n , on dn, (1.10) 

div(PV 2 w) -n+ ((PV 2 w)n-r), s = (M T ), S , onffi, (1.11) 

w = 0, on dD, (1.12) 

f = 0, on dD, (1.13) 

coupled with the equilibrium conditions (jl.6p . which has a unique solution 
w G H 2 (Q \ D). Therefore, the uniqueness question may be rephrased as 
follows: 

Given two solutions Wi to (I1.9p - (ll.l3p . ( II. 6p for D = D i} i = 1,2, satis- 
fying, for some g G A 

W!-w 2 = g, (w 1 -w 2 ), n = g, n , onT, (1.14) 

does D\ = D 2 hold? 

Obviously, the inverse problem above is equivalent to the determination 
of the portion dD of the boundary of Q \ D in the boundary value problem 
CT-flTT31). 

In the stability issue, instead of (11.71) and (jl.8p . we assume 

min|||wi -w 2 - g\\ L 2 (r) + \\( Wl - w 2 - g), n \\ L 2 (r) j < e, (1.15) 

for some e > 0, and we ask for the following estimate 

du(9Di,dIh)<v(e), (1.16) 
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where 77(e) is a suitable infinitesimal function. 

The uniqueness for the problem above has been proved in [M-R3j under 
the a priori assumption of C 3 ' 1 regularity of dD and with only one nontrivial 
couple field. Here, by nontrivial we mean that 



Concerning the stability issue, in |M-R-Ve5] we have proved a log-log 
type estimate, namely in inequality f l 1 . 1 6 j) we have r/(e) = O ((log | loge|) _Q ), 
where the positive parameter a depends on the a priori data, see Theorem 
14.31 below for a precise statement. 

In |M-R-Ve2] the inverse problem of determining a cavity in an elastic 
plate has been faced. We recall that in such conditions (1.4)-(1.5) are 

replaced by homogeneous Neumann boundary conditions, which are much 
more difficult to handle with respect to Dirichlet boundary conditions aris- 
ing in the case of rigid inclusions. For this reason a uniqueness result has 
been established making two linearly independent boundary measurements. 
In |M-R3j it has also been proved a uniqueness result for a variant of the 
problems considered above, that is the case of a plate whose boundary has 
an unknown and inaccessible portion where M = 0. In this case, thanks to 
the more favorable geometric situation, one measurement suffices to detect 
the unknown boundary portion. The corresponding stability results for these 
two cases have not yet been proved. 

The methods used to prove the above mentioned uniqueness and sta- 
bility results are based on unique continuation properties and quantitative 
estimates of unique continuation for solutions to the plate equation ( II. ip . 
Since such properties and estimates are consequences of the three sphere in- 
equality for solutions to equation (II. ip . we will discuss a while about the 
main features of such inequality. 

The three sphere inequality for solutions to partial differential equations 
and systems has a long and interesting history that intertwines with the issue 
of unique continuation properties and the issue of stability estimates |A1-M] . 
[Ho85] . [M], [J3], [LiE], |Lav-Rom-S] . |L-N-W] . |L-Nak-Wj . In many 

important cases, the three sphere inequality is the elementary tool to prove 
various types of quantitative estimates of unique continuation such as, for 
example, stability estimates for the Cauchy problem, smallness propagation 
estimates and quantitative evaluation of the vanishing rate of solutions to 
PDEs. Such questions have been intensively studied in the context of second 
order equations of elliptic and parabolic type. We refer to |Al-Ro-R-Vej and 
|Ve2j where these topics are widely investigated for such types of equations. 

The three sphere inequality for equation ( II .ip has been proved in |M-R-Ve5] 
under the very general assumption that the elastic material of the plate is 




(1.17) 
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anisotropic and obeys the so called dichotomy condition. Roughly speaking, 
such a condition implies that the plate operator at the left hand side of (II. ip 
can be written as L 2 L\ + Q, where L 2 , L\ are second order elliptic operators 
with C 1 ' 1 coefficients and Q is a third order operator with bounded coeffi- 
cients. For more details we refer to (I4.5ap - (14.5bl) below and |M-R- Ve5j . A 
simplified version of such inequality is the following one 

/ \V 2 w\ 2 <c(f \V 2 wA ( [ \V 2 wA , (1-18) 

JB T2 (x ) \jB ri {x ) J \JB rz {x ) J 

for every r\ < r 2 < r 3 , where 5 G (0, 1) and C depend only on the parameters 
related to the regularity, ellipticity and dichotomy conditions assumed on C, 
and on the ratios ri/r 2 , r 2 /r 3 ; in particular, 5 and C do not depend on w. 

Previously, the three sphere inequality was proved in |M-R-VeT] . for the 
isotropic plate (that is C m (x) = 5 ij 5 k iX(x) + (5 ik 5 j i + 5 u 5 jk ) /jl(x), i,j, k, I = 1, 2) 
and in |Gej for the class of fourth (and higher) order elliptic equation Cu = 
where C = L 2 L\ and L 2 , L\ are second order elliptic equation with C 1,1 
coefficients. 

The proof of the stability result, of which we give a sketch in Subsection 
14.21 has essentially the same structure of the proofs of analogous stability 
results in the following context: 

a) Second order elliptic equations: |Al-Roj . [Be- Vej . |Roj (two variables 
elliptic equations); |A1-B-R-Velj . [Al-B-R-Ve2] . [Si] (several variables elliptic 
equations) 

b) Second order parabolic equation: |C-R-Vel] . [C-R-Ve2] . [Dc-Ro-Ve] . 
[B-Dc-Si-Ve] . [VeT] . jVe2] 

c) Elliptic systems: [M-R3j . [M-R4j (elasticity); [BaJ (Stokes fluid) 

It is important to say that the stability estimates proved in the papers of 
list a) and b) are of logarithmic type, that is an optimal rate of convergence, 
as shown by counterexamples ( [All] . |Dc-Ro] for case a) and [Dc-Ro-Ve] . 
|Ve2j for case b)). The stability estimates proved in the papers of list c) 
and in the case of plate equation are of log-log type, that is with a worse 
rate of convergence. It seems difficult to improve such an estimate. We 
believe that the main difficulty to get such an improvement is due to the 
lack of quantitative estimates of strong unique continuation property at the 
boundary. In order to give an idea of the crucial point which marks the 
difference between these cases, let us notice that, by iterated application of 
the three sphere inequality (11.181) it can be proved there exists p > such 
that for every p G (0, p) and every x G dDj, j = 1, 2, the following inequality 
holds true 

/ _ \V 2 w 3 \ 2 > Cexp (~Ap- B ) , (1.19) 
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where A > 0, B > and C > only depend on the a priori information, in 
particular they depend by the quantity (frequency) 

\\M\\ L 2( d Q iR 2) 

\\M\\ i 

In cases a) and b), it is possible to prove a refined form of inequality (11.191) . 
in which the exponential term is replaced with a positive power of p, obtain- 
ing a quantitative estimate of strong unique continuation property at the 
boundary. It has been shown in [Al-B-R-Vel] that this is a key ingredient 
in proving that the stability estimate for the corresponding inverse problem 
with unknown boundaries in the conductivity context is not worse than log- 
arithm. This mathematical tool is available for second order elliptic, |A-E] . 
and parabolic equations, |Es-Fe-Ve] . but is not currently available for ellip- 
tic systems and plate equation. This happens even in the simplest case of 
isotropic material, and this is the reason for the presence of a double loga- 
rithm in our stability estimate. Finally, as remarked in |M-R-Ve4] . it seems 
hopeless the possibility that solutions to ( II. ip can satisfy even a strong unique 
continuation property in the interior, without any a priori assumption on the 
anisotropy of the material, see also [Ali] . Regarding this point, our dichotomy 
condition (I4.5a|) - fl4.5b[) basically contains the same assumptions under which 
the unique continuation property holds for a fourth order elliptic equation in 
two variables. 

In the present paper we have also proved constructive upper and lower 
estimates of the area of a rigid inclusion or of a cavity, D, in terms of an 
easily expressed quantity related to work. More precisely, suppose we make 
the following diagnostic test. We take a reference plate, i.e. a plate without 
inclusion or cavity, and we deform it by applying a couple field M at the 
boundary dfl. Let Wo be the work exerted in deforming the specimen. Next, 
we repeat the same experiment on a possibly defective plate. The exerted 
work generally changes and assumes, say, the value W. We are interested 
in finding constructive estimates, from above and from below, of the area 
of D in terms of the difference \W — Wq\. In order to prove such estimates 
we proceed along the path outlined in |M-R-VeT] and [M-R-Ve6] in which 
the inclusion inside the plate is made by different elastic material. In this 
introduction we illustrate such intermediate case, since the scheme of the 
mathematical procedure is fairly simple to describe. With regard to this 
intermediate case we also want to stress that, in contrast to the extreme 
cases, there are not available any kind of uniqueness result for the inverse 
problem of determining inclusion D from the knowledge of a finite number of 
measurements on the boundary . This appears to be an extremely difficult 
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problem. In fact, despite the wide research developed in this field, a general 
uniqueness result has not been obtained yet even in the simpler context 
arising in electrical impedance tomography (which involves a second order 
elliptic equation), see, for instance, |Isl] and |A11] for an extensive reference 
list. 

Denoting, as above, by w the transversal displacement of the plate and by 
M T , M n the twisting and bending components of the assigned couple field M, 
respectively, the infinitesimal deformation of the defective plate is governed 
by the fourth order Neumann boundary value problem 

~ div (div (( X n\ D W + Xd&)V 2 w)) = 0, inn, (1.21) 

(FV 2 w)n ■ n = -M n , on dQ, (1.22) 

dw(¥V 2 w)-n + ((¥V 2 w)n-T), s = (M r ), s , on dtt. (1-23) 

In the above equations, xd denotes the characteristic function of D. The 
plate tensors P, P are given by 

h 3 ~ h 3 ~ 
P=— C, P = — C, (1.24) 
12 ' 12 ' K J 

where C is the elasticity tensor describing the response of the material in the 
reference plate Q and satisfies the usual symmetry conditions (14. 2p . regularity 
condition (14.31) . strong convexity condition f)4.4p and the dichotomy condition, 
whereas C denotes the (unknown) corresponding tensor for the inclusion D. 
The work exerted by the couple field M has the expression 



W = - / M T , s w + M n w m . (1.25) 
Jon 

When the inclusion D is absent, the equilibrium problem fll.2ip - fll.23p be- 
comes 

div (div (PV 2 w )) = 0, in Q, (1.26) 

(PV 2 w )« • n = -M n , on <9ft, (1.27) 

div(PV 2 w ) -n+ ((PV 2 w )n-r)„= (M T ), S , onffi, (1-28) 

where Wq is the transversal displacement of the reference plate. The corre- 
sponding work exerted by M is given by 

W = - I M TiS w + M n w ,n. (1-29) 
Jan 
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In |M-R-Ve6] the following result has been proved. Assuming that the 
following fatness- condition is satisfied 



area ({a; e D\ dist{x,dD} > hi}) > iarea(D), (1.30) 



where hi is a given positive number, then 



Ci 



w -w. 



< area(L>) < C 2 



W-W ( 



;i.3i) 



where the constants Ci, C 2 only depend on the a priori data. Besides the as- 
sumptions on the plate tensor C given above, estimates (ll.3ip are established 
under some suitable assumption on the jump C — C. 

In the extreme cases corresponding to a rigid inclusion or a cavity D, 
the estimates are a little more involute than (ll.3ip and additional regularity 
conditions on the boundary 3D and on the plate tensor P are generally 
required. The main difference between extreme and intermediate cases lies 
in the estimate from below of |D|. Indeed, in the former case we use regularity 
estimates (in the interior) for the reference solution wq to equation (11.261) . 
whereas in the latter we combine regularity estimates, trace and Poincare 
inequalities. On the other hand, the argument used for the estimate from 
above of \D\ is essentially the same as in |M-R-Ve6] and involves quantitative 
estimates of unique continuation in the form of three sphere inequality for 
the hessian V 2 u>o- It is exactly at this point that the dichotomy condition 
(I4.5ap - fl4.5bl) on the tensor C is needed. 

The analogous bounds in plate theory for intermediate inclusions were 
first obtained when the reference plate satisfies isotropic conditions, |M-R- VeT] . 
extended to anisotropic materials satisfying the dichotomy conditions, |M-R-Ve6] , 
obtained in a weaker form for general inclusions in absence of the fatness con- 
dition, |M-R-Ve3] , and recently in the context of shallow shells in |Dc-Li- Wa] 
and jDc-Li-Ve-Wa] . The reader is referred to |K-S-Shj . [AT-R] . |A1-R-Sj . 
|Be-Fr-Ve] for size estimates of inclusions in the context of the electrical 
impedance tomography and to [Ik], |A1-M-R2] . |A1-M-R3] . LM-M-R-V for 
corresponding problems in two and three-dimensional linear elasticity. See 
also |L-Dc-Nj for an application of the size estimates approach in thermog- 
raphy. 

Size estimates for extreme inclusions were obtained in |A1-M-Ri] for elec- 
tric conductors and in |M-R1] for elastic bodies, see also |A1-M-R2] . 

The paper is organized as follows. In Section [2] we collect some notation. 
In Section [3] we provide a derivation of the Kirchhoff-Love model of the plate. 
In Section @] we present the uniqueness and stability results concerning the 
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determination of rigid inclusions, and the uniqueness results for the case of 
cavities and unknown boundary portions. In particular we have focused our 
attention on the case of rigid inclusions and, for a better comprehension of 
the arguments, we have recalled the proof of the uniqueness result, using it as 
a base for a sketch of the more complex proof of the stability result. Section 
contains the estimates of the area of extreme inclusions. 

2 Notation 

Let P = (x 1 (P),x 2 (P)) be a point of R 2 . We shall denote by B r (P) the disk 
in R 2 of radius r and center P and by R a ^{P) the rectangle R a ^(P) = {x = 
(xi,x 2 ) | \x\ — xx(P)\ < a, \x 2 — x 2 (P)\ < b}. To simplify the notation, we 
shall denote B r = B r (0), R a ,b = R a ,b{0). 

Definition 2.1. (C jl regularity) Let Q be a bounded domain in R 2 . Given 
k G N, we say that a portion S of dQ is of class C ,x with constants po, 
Mo > 0, if, for any P G S, there exists a rigid transformation of coordinates 
under which we have P = and 

^ n#j£_ - = {x = (xx,x 2 ) G Rpo p | x 2 >^(a;i)}, 

where ij) is a C k ' 1 function on ^— j^, -^J satisfying 

^(0) = 0, ip'(0) = 0, when k > 1, 

When k = we also say that S is of Lipschitz class with constants po, Mq. 

Remark 2.2. We use the convention to normalize all norms in such a way that 
their terms are dimensionally homogeneous with the L°° norm and coincide 
with the standard definition when the dimensional parameter equals one, see 
|M-R-Ve5j for details. 

Given a bounded domain Q in R 2 such that dQ is of class C ,1 , with 
k > 1, we consider as positive the orientation of the boundary induced by 
the outer unit normal n in the following sense. Given a point P G dfl, let 
us denote by r = r(P) the unit tangent at the boundary in P obtained by 
applying to n a counterclockwise rotation of angle ?, that is r = e% x n, 
where x denotes the vector product in R 3 , {ei,e 2 } is the canonical basis in 
R 2 and e% = e\ x e 2 . Given any connected component C of dQ and fixed 
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a point P G C, let us define as positive the orientation of C associated to 
an arclength parametrization (p(s) = (x%(s), x 2 (s)), s G [0, 1(C)], such that 
(p(0) = P and ip'(s) = r(ip(s)), where 1(C) denotes the length of C. 

Throughout the paper, we denote by diii, d s u, and d n u the derivatives of 
a function u with respect to the Xi variable, to the arclength s and to the 
normal direction n, respectively, and similarly for higher order derivatives. 

We denote by M 2 the space of 2 x 2 real valued matrices and by C(X, Y) 
the space of bounded linear operators between Banach spaces X and Y. 

For every 2x2 matrices A, B and for every L G £(M 2 , M 2 ), we use the 
following notation: 

(1^4)^' = Lij k iA k i, (2.1) 
A - B = AijBij, \A\ = (A-A) 1 i. (2.2) 

Notice that here and in the sequel summation over repeated indexes is im- 
plied. 

Finally, let us introduce the linear space of the affine functions on R 2 
A = {g(x\, X2) = ax\ + bx 2 + c, a, b, c G M}. 



3 The Kir chhoff- Love plate model 

In the last two decades different methods were used to provide new justifi- 
cation of the theory of thin plates. Among these, we recall the method of 
asymptotic expansion |C-Dj . the method of internal constraints [PG| . [L-PGj . 
the theory of T-convergence in conjunction with appropriate averages |A-B-Pj 
or on rescaled domain and with rescaled displacements |B-C-G-R] , |Pal] , and 



weak convergence methods on a rescaled domain and with rescaled displace- 
ments [C] . We refer the interested reader to |Pa2] for a recent account of the 
advanced results on this topic. The present section has a more modest aim: 
to show how to deduce the equations governing the statical equilibrium of an 
elastic thin plate following the classical approach of the Theory of Structures. 

Let us consider a thin plate f2 x [— | , |] with middle surface represented 
by a bounded domain Q in M 2 having uniform thickness h, h <<diam(fi), 
and boundary dfl of class C 1 ' 1 . Only in this section, we adopt the convention 
that Greek indexes assume the values 1,2, whereas Latin indexes run from 1 
to 3. 

We follow the direct approach to define the infinitesimal deformation of 
the plate. In particular, we restrict ourselves to the case in which the points 
x = (xi,x 2 ) of the middle surface Q are subject to transversal displace- 
ment w(xi, x 2 )e 3 , and any transversal material fiber {x} x [— |, |], x G ft, 
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undergoes an infinitesimal rigid rotation oo(x), with u(x) ■ e 3 — 0. In this 
section we shall be concerned exclusively with regular functions on their do- 
main of definition. For example, the above functions w and u are such that 
w G C°°(f2,R) and u> G C°°(fi,R 3 ). These conditions are unnecessarily re- 
strictive, but this choice simplifies the mechanical formulation of the equilib- 
rium problem. The above kinematical assumptions imply that the displace- 
ment field present in the plate is given by the following three-dimensional 
vector field: 

u(x,x 3 ) = w(x)e 3 + x 3 ip(x), x G f2, \x 3 \ < ^, (3.1) 

where 

<f(x) = u(x) x e 3 , x G Q. (3.2) 

By (13. ip and (13. 2p . the associated infinitesimal strain tensor E[u] G M 3 takes 
the form 

E[u](x,x 3 ) = (Vn)^ m (x,x 3 ) = x 3 (VMx)) sym + ® e 3 ) sym , (3.3) 

where V x (-) = af~(') e « * s the surface gradient operator, V SJ/m (-) = |(V(-) + 
V T (-)), and 

l(x) = ip(x) + V x w(x). (3.4) 

Within the approximation of the theory of infinitesimal deformations, 7 is the 
angular deviation between the transversal material fiber at x and the normal 
direction to the deformed middle surface of the plate at x. In Kirchhoff-Love 
theory it is assumed that every transversal material fiber remains normal to 
the deformed middle surface, e.g. 7 = in Q. 

The traditional deduction of the mechanical model of a thin plate follows 
essentially from integration over the thickness of the corresponding three- 
dimensional quantities. In particular, taking advantage of the infinitesimal 
deformation assumption, we can refer the independent variables to the initial 
undeformed configuration of the plate. 

Let us introduce an arbitrary portion tt'x [— |, |1 of plate, where CC 
is a subdomain of Q with regular boundary. Consider the material fiber 
{x} x [— |, |] for x G dQ' and denote by t(x,x 3 ,e a ) G M 3 , \x 3 \ < f, the 
traction vector acting on a plane containing the direction of the fiber and 
orthogonal to the direction e a . By Cauchy's Lemma [Tj we have t(x, x 3 , e a ) = 
T(x,x 3 )e a , where T(x,x 3 ) G M 3 is the (symmetric) Cauchy stress tensor at 
the point (x,x 3 ). Denote by n the unit outer normal vector to dQ' such 
that n ■ e 3 = 0. To simplify the notation, it is convenient to consider n as a 
two-dimensional vector belonging to the plane x 3 = containing the middle 
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surface Q of the plate. By the classical Stress Principle for plates [VI] . we 
postulate that the two complementary parts Q' and fi \ fi' interact with one 
another through a filed of force vectors R = R(x, n) G M 3 and couple vectors 
M = M(x,n) G M 3 assigned per unit length at x G dQ'. Denoting by 

rh/2 

R{x 1 e a ) = I t(x, X3, e a )dx3 (3.5) 

J-h/2 

the force vector (per unit length) acting on a direction orthogonal to e a and 
passing through x G dQ', the contact force -R(x, n) can be expressed as 

R(x,n) = T n (x)n, x G dtt', (3.6) 

where the surface force tensor T Q (x) G M 3x2 is given by 

T n (x) = R(x, e a ) g> e Q , in fl (3.7) 

Let P = / — e% g) e3 be the projection of M 3 along the direction e^. T Q is 
decomposed additively by P in its membranal and shearing component 

r-pQ prpQ _|_ p^rpQ — rpQ(m) _|_ rpQ(s) 

where, following the standard nomenclature in plate theory, the components 
^ap (= ), a,(3 = 1,2, are called the membrane forces and the com- 
ponents T^ s \ (3 = 1,2, are the shear forces (also denoted as Tfg = Qp)- 
The assumption of infinitesimal deformations and the hypothesis of vanishing 
in-plane displacements of the middle surface of the plate allow us to take 

T n{m) = 0> in Q ^ 

Denote by 



rh/2 

M{x 1 e a ) = I x-jCz x t(x, X3, e a )dx3, a = 1,2, (3.10) 

J-h/2 



the contact couple acting at x G dVt' on a direction orthogonal to e Q passing 
through x. Note that M(x, e a )-e^ = by definition, that is M(x, e a ) actually 
is a two-dimensional couple field belonging to the middle plane of the plate. 
Analogously to (13.61) . we have 

M(x,n) = M n (x)n, x G dQ' , (3.11) 

where the surface couple tensor M n (x) G M 2x2 has the expression 

M n (x) = M(x,e a ) ®e a . (3.12) 
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A direct calculation shows that 



M(x,e a ) = e 3 x epMp a (x), 



(3.13) 



where 



M Pa (x) 



h/2 



x 3 Tp a (x, x 3 )dx 3 , a, (3 = 1, 2, 



(3.14) 



-h/2 



are the bending moments (for a = 0) and the twisting moments (for a 7^ /3) 
of the plate at x (per unit length). 

The differential equilibrium equation for the plate follows from the inte- 
gral mechanical balance equations applied to any subdomain Q' CC Q [T]. 
Denote by q(x)e 3 the external transversal force per unit area acting in Q. 
The statical equilibrium of the plate is satisfied if and only if the following 
two equations are simultaneously satisfied: 

Ian' T n nds + j n , qe 3 dx = 0, (3.15) 
Idn' (( x — x o) x + M n nj ds + J n ,{x — x ) x qe 3 dx = 0, (3.16) 

for every subdomain Q' C Q, where xq is a fixed point. By applying the 
Divergence Theorem in Q,' and by the arbitrariness of Q' we deduce 



div^T^) + qe 3 = 0, in Q, 

div x M n + (T n ^) T e 3 x e 3 = 0, in fi. 



(3.17) 
(3.18) 



Consider the case in which the boundary of the plate dfl is subjected si- 
multaneously to a couple field M, M ■ e 3 = 0, and a transversal force field 
Qe 3 . Local equilibrium considerations on points of dfl yield the following 
boundary conditions: 



M n n = M, on dtt, 
T n ^n = Qe 3 , onffi. 



(3.19) 
(3.20) 



where n is the unit outer normal to dQ. In cartesian components, the equi- 
librium equations (I3.17p -f l3.20p take the form 



M a p t /3 — Q a — 0, 

Qa,a + q = 0, 

M a/3 n a np = M n , 
M a/3 T a np = -M T , 



in Q, a = 1, 2, 
in Q, 
on dfl, 
on dfl, 
on <9fl 



(3.21) 
(3.22) 
(3.23) 
(3.24) 
(3.25) 
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Here, following a standard convention in the theory of plates^jve have decom- 
posed the boundary couple field M in local coordinates as M = M T n + M n r. 

To complete the formulation of the equilibrium problem, we need to in- 
troduce the constitutive equation of the material. We limit ourselves to the 
Kirchhoff-Love theory and we choose to regard the kinematical assumptions 
-^i3[ M ] — 0, i = 1,2,3 (see (13. 3p . with 7 = 0) as internal constraints, that 
is we restrict the possible deformations of the points of the plate to those 
whose infinitesimal strain tensor belongs to the set 

M = {E G M 3x3 \E = E T , E ■ A = 0, for A = e { <g> e 3 + e 3 <g> e h i = 1,2, 3}. 

(3.26) 

Then, by the Generalized Principle of Determinism [Tj, the Cauchy stress 
tensor T at any point (x, x 3 ) of the plate is additively decomposed in an 
active (symmetric) part T A and in a reactive (symmetric) part Tr: 

T = T A + T R , (3.27) 

where Tr does not work in any admissible motion, e.g., Tr G M. . Consis- 
tently with the Principle, the active stress T A belongs to Ai and, in cartesian 
coordinates, we have 

T A = T Aa pe a ®e p , a, (3 = 1,2, (3.28) 

T R = T Ra3 e a <g> e 3 + T ma e 3 ® e a + T m3 e 3 g> e 3 . (3.29) 

In linear theory, on assuming the reference configuration unstressed, the ac- 
tive stress in a point (x,x 3 ) of the plate, x G Q and |x 3 | < h/2, is given by 
a linear mapping from Ai into itself by means of the fourth order elasticity 
tensor Cj^: 

T A = C M E[u}. (3.30) 

Here, in view of fl3T2"6D and (13381 . C M can be assumed to belong to £(M 2 , M 2 ). 
Moreover, we assume that Cm is constant over the thickness of the plate and 
satisfies the minor and major symmetry conditions expressed in cartesian 
coordinates as (we drop the subscript Ai) 

Cat3~j5 = Cpa-yS = Ca/3Sf = C-ySa/3, OC, /3,J,5 — 1,2, in Q. (3.31) 

We refer to |PG] and |L-PGj for a representation formula of C based on the 
maximal response symmetry of the material compatible with the internal 
constraints. 

Using (13.271) and recalling (13.91) . we obtain the corresponding decomposi- 
tion for T n and M n : 

T n = T^ {s \ M n = M%, (3.32) 
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that is the shear forces and the moments have reactive and active nature, re- 
spectively. By (13.301) . after integration over the thickness, the surface couple 
tensor is given by 

U3 

M n {x) = - — £C(x)(V>(x)), in tt, (3.33) 

where £ G M 2 has cartesian components E\\ = E\\ = 0, £u = —1, £21 = 1- 
Constitutive equation (13.331) can be written in more expressive way in terms 
of the bending and twisting moments as follows: 

M a p(w) = -P a ^ s (x)w, lS , «, = 1,2, (3.34) 

where 

F(x) = y^C(x), in Q, (3.35) 

is the plate elasticity tensor. Combining (I3.17P and (I3.18p . and by eliminating 
the reactive term T^ s \ we obtain the classical partial differential equation 
of the Kirchhoff- Love's bending theory of thin elastic plates, that, written in 
cartesian coordinates, takes the form 

{P a /3js(x)w n s {x)) ia p = Q, in ^- (3-36) 

In the remaining part of this section we complete the formulation of the 
equilibrium problem for a Kirchhoff-Love plate by writing the boundary con- 
ditions corresponding to ( 13. 230 - 03.250 . The determination of these bound- 
ary conditions is not a trivial issue because, first, a constitutive equation for 
shear forces is not available since these have a reactive nature (see (I3.32p ). 
and, second, because the three mechanical boundary conditions f l3.23j> — fl3.25j) 
should reasonably "collapse" into two independent boundary conditions for 
the fourth order partial differential equation (13.361) . To this aim, under the 
additional assumption of C positive definite, we adopt a variational approach 
and we impose the stationarity condition on the total potential energy func- 
tional J of the plate. Consider the space of regular kinematically admissible 
displacements 

V = {v: Qx(-h/2, h/2) M 3 | v(x, x 3 ) = r](x)e 3 -x 3 V x r](x), with rj : K}. 

(3.37) 

The energy functional J : T> — V R is defined as 

J(v) = a{v,v) -l{v), (3.38) 

where a(v,v) is interpreted as the elastic energy stored in the plate for the 
displacement field v and l(v) is the load potential that accounts for the energy 
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of the system of applied loads q, M, Q. The three-dimensional expression of 
the elastic energy in the Linear Theory of Elasticity is given by 



h/2 

h/2 Jn 



a(v,v) = i / / CE[r] ■ E[r}<!.r<1.r :i . (3. 39) 



where, in view of (13.371) . E[v] = — x 3 V^(x). After integration over the 
thickness, we obtain 



a(v,v) = -- I M a p(r])ri, a pdx, (3.40) 



where M a p(rf) are as in (I3.34[) with w replaced by rj. The load functional has 
the expression 

l(v) = / q7]dx + {Qr] + M 2 ri,2 -Mi^,i )ds. (3.41) 
Jn Jan 

Then, the stationarity condition (in fact, minimum condition) on J at w 
yields 

/ M a p(w) V , a pdx+ I qr)+ {Q7 ] + M 2 7 ] , 2 -M 1 r ] , l )ds = 0, (3.42) 
Jn Jn Jan 

for every regular function rj. Integrating by parts twice on the first integral 
we obtain 

/ (M aPy0l p(w) + q)r]dx + / {-M aj3 fi(w)n a + Q)t]ds+ 
Jn Jan 



+ / (M a p(w)npr 1a + M 2 T], 2 -M 1 ri, 1 )ds = 0. (3.43) 
Jan 

We elaborate the last integral I\ of (I3.43[) by rewriting the first order deriva- 
tives of rj on dQ in terms of the normal and arc-length derivative of 77. We 
have 

h= (M a i3(w)n/3n a - M x n x + M 2 n 2 )?7, n ds+ 
Jan 

+ / (M a p(w)n /3 T a + M 1 n 2 + M 2 n 1 )ri, s ds = + (3.44) 
Jan 

and, integrating by parts on dfl, we get 
I[ = (M a p(w)npr a + M x n 2 + M 2ni ) V \ s = l Q {dn) - 

{M aP {w)npT a + M x n 2 + M 2 m), s r]ds, (3.45) 



an 



16 



where l(dQ) is the length of dfl. Since dfl is of class C 1 ' 1 , the boundary 
term on the right end side of (I3.45j) identically vanishes. Therefore, the 
stationarity condition of J at w takes the final form 




/ [-(M a p(w)npT a ), s -M a p t p(w)n a + Q- {M x n 2 + M 2 n 1 ), s ) r]ds+ 
an v J 



(M a/3;a/3 (w) + q)rjdx+ 




Jon 



for every rj G C°°(fi, R). By the arbitrariness of the function rj, and of the 
traces of r] and r], n on dfl, we determine the equilibrium equation (I3.36P and 
the desired Neumann boundary conditions on dQ: 



4 Uniqueness and stability of extreme inclu- 
sions and free boundaries 

4.1 Rigid inclusions: uniqueness 

In the sequel we shall assume that the plate is made of nonhomogeneous 
linear elastic material with plate tensor 

P=^C, (4.1) 

where the elasticity tensor C(x) G £(M 2 , M 2 ) has cartesian components Cijki 
which satisfy the following symmetry conditions 



M aj3 {w)n a np = M, 



(3.47) 



M a p )P {w)n a + (M a p(w)npT a ), a = Q- (M T ) 



(3.48) 



C^ki — Ckuj — Ckiji k, I — 1, 2, a.e. in Q, 



(4.2) 



and, for simplicity, is defined in all of M. 2 . 
We make the following assumptions: 
I) Regularity 



C G C7 1,1 (R 2 ,£(M 2 ,M 2 )) J 



(4.3) 



II) Ellipticity (strong convexity) There exists 7 > such that 
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CA-A > 7 |A| 2 , 

for every 2x2 symmetric matrix A. 
Ill) Dichotomy condition 



in 



either V(x) > 0, for every x G 
or T>(x) = 0, for every x G 



where 



a 







V 


(*) = 


1 | 

a 


det£( 


x)|, 








/ 


a 


Ai 


«2 




04 





\ 









a 




«2 




a,{ 















a 


a l 


0-2 


a 3 


04 


S(x) = 




4a 


3ai 


2a 2 


a 3 


















4a 


3ai 


2a 2 


03 


















4a 


3ai 


2a 2 


a 3 







V 











4a 


3ai 


2a 2 


a-3 J 


A , ai = 




4C , 


a 2 = 


2B + 4£ , a 3 


= 4L> 


0, 04 



and 



C1111 


= 4>, 


C1122 = 


C2211 = Bo 


• 


C1112 


= C1121 


= C1211 


— C2111 = 


Co 


< C2212 


= C2221 


= C1222 


= C2122 = 


£0 


C1212 


= C1221 


= C2112 


= C2121 


E 


C2222 


= F , 









Remark 4.1. Whenever (I4.5ap holds we denote 

5\ = minD. 



(4.4) 

(4.5a) 
(4.5b) 

(4.6) 



(4.7) 



(4.i 



(4.9) 



(4.10) 



We emphasize that, in all the following statements, whenever a constant is 
said to depend on 5± (among other quantities) it is understood that such 
dependence occurs only when f)4.5al) holds. 

On the assigned couple field M let us require the following assumptions: 
M G H-i(dn,R 2 ), (M n ,M r , s )^0, (4.11) 



Mi = 0, 



1,2. 



(4.12) 



<■■><> 
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Theorem 4.2 (Unique determination of a rigid inclusion with one measure- 
ment). Let Q be a simply connected domain in R 2 such that dQ is of class C 1 ' 1 
and let D i; i = 1,2, be two simply connected domains compactly contained 
in Q, such that dDi is of class C 3,1 , i = 1,2. Moreover, let T be a nonempty 
open portion of dQ, of class C 3,1 . Let the plate tensor P be given by (14. ip . 
and satisfying (14.21) - (14.41) and the dichotomy condition (I4.5ap or (14.5b[) . Let 
M be a boundary couple field satisfying (14.111) - (I4.12p . Let Wi, i = 1, 2, be the 
solutions to the mixed problem (ll.9p - (11.131) . coupled with (II ,6p . with D = Di. 
If there exists g G A such that 

w 1 -w 2 = g, {w 1 -w 2 ), n = g,n, onT, (4.13) 

then 

D 1 = D 2 . (4.14) 

Proof of Theorem \4-S\ Let G be the connected component of Q \ (D\ U D 2 ) 
such that T C dG. Let us notice that, since w,i satisfies homogeneous Dirich- 
let conditions on the C 3,1 boundary dDi, by regularity results we have that 
Wi G H 4 (Q \ Di), for every Q, Di CC Q CC Q, i — 1,2 (see, for example, 
|Ag|). By Sobolev embedding theorems (see, for instance, [Ad]), we have 
that Wi and Vwi are continuous up to dDi, % = 1,2. Therefore 

u>i = 0, V< = 0, on^A- (4.15) 

Let w = wi — w 2 — g, with g{x\, x 2 ) = ax\ + bx 2 + c. By our assumptions, w 
takes homogeneous Cauchy data on T. From the uniqueness of the solution to 
the Cauchy problem (see, for instance, Theorem 3.8 in |M-R-Ve4] ) and also 
Remark 4 in |M-R-Ve2] ) and from the weak unique continuation property, 
we have that 

w = 0, in G. 
Let us prove, for instance, that D 2 C D\. We have 

D 2 \D~ 1 cn\(D 1 UG), 



d{Q \ (D 1 U GO) = Ei U E 2 , 

where S 2 = dD 2 n dG, E x = d(Q \ (AUG)) \ S 2 C dD x . 
Let us distinguish two cases 

i) dD 1 n E 2 ± 0, 

ii) dD l n E 2 = 0. 

In case i), there exists Po £ <9-A H S 2 . Since Wi(P ) = 0, w(Po) = 0, we 
have g(Po) = 0. 
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Let P n 6 G, P n -)> P Q . We have 

Vw(P n ) = 0, 

= lim Vw(P n ) = V<(P ) - V^(Pq) - (a, 6) = -(a, 6), 

n— >oo 

so that g = c, and, since g(Po) = 0, 

. „ f u>i = w 2 = 0, on £ 2 , 

g = ffliEto 2 , mG, => <^ _ (4.16) 



Vu>i = Vu? 2 = 0, on £ 2 . 
Integrating by parts equation (div(div(PV 2 wi)))wi = 0, we obtain 



/ 

J a 



¥V 2 Wi ■ W = / (PV^wiv ■ u)w 1;n + 



n\(r>iuG) JEius 2 

+ / (PVV^r)w M - / (div(PVV) ■ (4.17) 



where z/ is the outer unit normal to f2 \ (Di U G). By 

U7i = 0, Vwi = 0, on Ei, 

and by f)4.16p . we have 

Wi = w 2 = 0, Vifi = Vw 2 = 0, on £ 2 , 

so that 

0=/ PVV-VV>7 / |VV| 2 - (4.18) 



./f2\CDiUG) JD 2 \D 1 

If D 2 \D 1 7^ 0, then, by the weak unique continuation principle, coincides 
with an affine function in Q\ Di, contradicting the choice of the nontrivial 
Neumann data M on dQ. Therefore D 2 C Di and, by the regularity of Di, 
1 = 1,2, D 2 cDx. _____ 

In case ii), either D x fl D 2 = or Di C -D 2 . 

Let us consider for instance the first case, the proof of the second case 
being similar. Integrating by parts, we have 

/ PVV ■ VV = / PVV ■ V 2 (wi -g) = 
Jd 2 Jd 2 

{FV 2 w 1 v-v){w 1 -g) m + 

dD 2 

I (pvV-t)^-^,-/ (div(PvV)-^K-^). 

>dD 2 JdD 2 



20 



By the regularity of dD 2 , we may rewrite it as 



/ PVV-VV= / (FV 2 w 1 is-is)(w 1 -g), n + 

JD 2 JdD 2 

- [ ((PVV*/ • t), 8 +div(PVV) • v) ( Wl - g), 
J a Do 



Recalling that wi — g = w 2 = 0, {w\ — g), n — W2,« = on dD 2 , we have 
that 

/ pv V • vV = o. 

Jd 2 

If D 2 7^ 0, then W\ coincides with an affine function in D 2 , and, by the 
weak unique continuation principle, also in Q\Di, contradicting the choice 
of a nontrivial M. Therefore D 2 = 0. Symmetrically, we obtain that D\ = 0, 
that is D 1 — D 2 . 

□ 

4.2 Rigid inclusions: stability 

In order to prove the stability estimates, we need the following further quan- 
titative assumptions. 

Given p , M , M 1 > 0, we assume that 

\n\ < M lP 2 , (4.19) 

dist(D,dn) > po, (4.20) 

<9f2 is of class C 2 ' 1 with constants po,M , (4-21) 

r is of class C 3 ' 1 with constants po, Mq, (4.22) 

dD is of class C 3 ' 1 with constants p , Mq, (4.23) 

where denotes the Lebesgue measure of Q. Moreover, we assume that for 
some P G £ and some S , < 5 < 1, 



^n^ (Po)cr, (4.24) 

and that 



|r| < (l-8 )\dn\. (4.25) 
On the Neumann data M we assume that 

Ma 2 (ffl,l 2 ), (M n ,(M T ), 8 ) ^0, (4.26) 
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Mi = 0, i = 1,2, (4.27) 

an 

supp(M) CC T, (4.28) 
and that, for a given constant F > 0, 



On the elasticity tensor C, we assume the same a priori information made 
in Subsection 14.11 and we introduce a parameter M > such that 

2 2 

Z)Po , HV m ai«IU-(R>)<M (4.30) 

i,j,k,l=l m=0 

We shall refer to the set of constants M , Mi, <5 , -F, 7, M, 5i as the a 
priori data. The scale parameter po w iH appear explicitly in all formulas, 
whereas the dependence on the thickness parameter h will be omitted. 

Theorem 4.3 (Stability result). Let Q be a bounded domain in M. 2 satisfying 
(I4.19p and (I4.2ip . Let Di, i — 1,2, be two simply connected open subsets of 
Q satisfying (I4.20p and (I4.23P . Moreover, let V be an open portion of dfl 
satisfying ( Q2J) , f04"P and fl425|) . Let M G L 2 (dn,R 2 ) satisfy ( 14^26]) - 
fT4T29]) and let the plate tensor P given by <KT$_ sa,tisfv (@~2]) ; fTOnjl . (Oft 
and the dichotomy condition. Let Wi G H 2 (Q \ Di) be the solution to (\1.9\i - 
fll,13p . coupled with ( II. 6p . u>/ien D = Di, i = 1,2. If, given e > 0, we 
have 

mm - W 2 -g\\ L \T,) +P0 IK^i ~W 2 -flO.nlljr^E)} ^ e > ( 4 - 31 ) 

then we have 

d H (dD l ,dD 2 ) <Cp Q (log|logel)^, 0<?<e-\ (4.32) 

and 

d H (D\,LT 2 ) < Cp (log|log?|)-", <?< e"\ (4.33) 
where e = 0u -+ € u and C, n, C > 0, < 77 < 1, are constants only 

pl\\M\\ _i ' " ' y 

depending on the a priori data. 
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Proof of Theorem \4-3\ Let us rough-out a sketch of the proof, referring the 
interested reader to Section 3 of [M-R- Ve5] . Retracing the proof of the 
uniqueness theorem, the basic idea is that of deriving the quantitative version 
in the stability context of the vanishing of the integral Jq w d uG \ PV 2 w« • V 2 u>«, 
for i = 1,2, that is a control with some small term emerging from the bound 
f!4.3ip on the Cauchy data. To this aim, we need stability estimates of con- 
tinuation from Cauchy data and propagation of smallness estimates. In par- 
ticular, the propagation of smallness of |V 2 w,|, for i — 1, 2, from a neighboor- 
hood of T towards d(Q \ (Di UG)) is performed through iterated application 
of the three sphere inequality fll . 18 j) over suitable chains of disks. A strong 
hindrance which occurs in this step is related to the difficulty of getting ar- 
bitrarily closer to the boundary of the set Q \ (D^ U G), due to the absence, 
in our general setting, of any a priori information on the reciprocal position 
of D\ and D 2 . For this reason, as a preparatory step, we derive the following 
rough estimate 



max 



2„,, |2 



tJD 2 \D 1 </DiYD 2 

which holds for every e < e~ x , with e 



VW^HMH^Qogllog^-s (4.34) 



and where C > depends 



only on 7, M, Si, M , Mi and S . 

Since a pointwise lower bound for |V 2 Wj| 2 cannot hold in general, the 
next crucial step consists in the following Claim. 
Claim. If 



max 



|vV| 2 , 



£>a\£>i 



£>i\£> 2 



V 



2 ■ 



(4.35) 



then 



d H (dD 1 ,dD 2 )<C Po 



log 



u " "i?-2(an,iR™) 
V 



(4.36) 



where B > and C > only depend on 7, M, Si, M , M 1 , S and F. 

Proof of the Claim. Denoting for simplicity d = du(dDi,dD-i), we may as- 
sume, with no loss of generality, that there exists xq G dD\ such that 
dist(x , dD 2 ) = d. Let us distinguish two cases: 

i) B d (x ) C D 2 ; 
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ii) B d (x ) n D 2 = 0. 

In case i), by the regularity assumptions made on dDi, there exists X\ G 
D 2 \ Di such that B td (xi) C D 2 \D X , with t = / 

1 +V 1+M o 

In |M-R-Ve6] . by iterated application of the three sphere inequality (ll.lSp . 
we have obtained the following Lipschitz propagation of smallness estimate: 
there exists s > 1, only depending on 7, M, Si, M and So, such that for 
every p > and every x £ D) sp , we have 



Bp{s) exp 



W > ,~ ru ^ ||M||^, (4.37) 



.4 



where A > 0, 5 > and C > only depend on 7, M, <5 X , M , Mi, S and F. 
By ( I4.35P and by applying (I4.37P with p = —, we have 

v> 4 (4.38) 



cxp 



where A > 0, 5 > and C > only depend on 7, M, Si, M , M x , S and F. 
By (Qg|| we easily find (OS]l . 

Case ii) can be treated similarly by substituting with w 2 - D 

By applying the Claim to f )4.34p . we obtain a first stability estimates of 
log-log-log type. At this stage, a tool which turns out to be very useful is 
a geometrical result, firstly stated in |A1-B-R-Vel] . which ensures that there 
exists e > 0, only depending on 7, M, Si, M , M 1; S and F, such that 
if e < e then dG is of Lipschitz class with constants po, L, with L and 
^ only depending on M . Lipschitz regularity prevents the occurrence of 
uncontrollable narrowings or cuspidal points in G and allows to refine the 
geometrical constructions of the chains of disks to which we apply the three 
sphere inequality, obtaining the better estimate 

max! / _|W| 2 , / _\V 2 w 2 \ 2 ) < p 2 Q \\M\\ 2 i|log?T CT , (4-39) 

which holds for every? < 1, where C > and o > depend only on 7, M, Si, 
Mo, Mi, So, L and ^. Again, by applying the Claim, the desired estimates 
follow. □ 



4.3 Cavities and unknown boundary portions: unique- 
ness 

In this subsection we consider the inverse problems of determining unknown 
boundaries of the following two kinds: i) the boundary of a cavity, ii) an 



24 



unknown boundary portion of dQ. In both cases, we have homogeneous 
Neumann conditions on the unknown boundary. Neumann boundary con- 
ditions lead to further complications in the arguments involving integration 
by parts. To give an idea of the differences, given any connected component 
F of D 2 \ Di, whose boundary is made of two arcs r C dD\ and 7 C dD 2 , 
having common endpoints Pi and P 2 , the analogue of (I4.17P becomes 

/ PVV • VV = / (PVVn 1 • rV), s - / (PV 2 w 2 n 2 • t 2 w 2 ), s = 

JF Jt J-y 

= [(PVVn 1 • r^iPx) - (FV 2 Wl n 2 ■ r 2 )(Px)] (w 1 (P 1 ) - w 1 {P 2 )), (4.40) 

where n l , r % denotes the unit normal and tangent vector to Di, i = 1,2. 

Since, in general, the boundaries of D\ and D 2 intersect nontangentially, 
the above expression does not vanish and the contradiction arguments fails. 
For this reason, we need two boundary measurements to prove uniqueness, 
as stated in Theorem 14.41 Instead, uniqueness with one measurement can 
be restored in the problem of the determination of an unknown boundary 
portion, by taking advantage of the fact that the two plates have a common 
regular boundary portion, say T, see Theorem 14.51 

Let D CC O be a domain of class C 1 ' 1 representing an unknown cavity 
inside the plate fl Under the same assumptions made in Subsection 14. 14 the 
transversal displacement w satisfies the following Neumann problem 



div(div(PV 2 w)) = 0, 


in n\D, 


(4.41) 


(PV 2 w)n ■ n = -M n , 


on dQ, 


(4.42) 


div(PV 2 w) ■ n + {(FV 2 w)n ■ t), s = (M T ), S , 


on dQ, 


(4.43) 


(PV 2 to)n • n = 0, 


on dD, 


(4.44) 


div(PV 2 w) • n + ((FV 2 w)n ■ t), s = 0, 


on dD, 


(4.45) 



V 



which admits a solution w G H 2 (Q \ Q), which is uniquely determined up to 
addition of an affine function. 

Concerning the inverse problem of the determination of the cavity D 
inside the plate, let us recall the following result. 

Theorem 4.4 (Uniqueness with two boundary measurements). Let Q be a 
simply connected domain in R 2 such that dQ is of class C 1,1 and let Di, 
i = 1,2, be two simply connected domains compactly contained in Q, such 
that dDi is of class C ,l , i = 1,2. Moreover, let T be a nonempty open portion 
of dVt, of class C 3 ' 1 . Let the plate tensor ^P be given by f !4.ip . and satisfying 
(14. 3p f)4.4p and the dichotomy condition (14 .5 aft or (jjJabp . Let M, M* be two 
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boundary couple fields both satisfying fl4.lip - fl4.12p and such that (M n ,M TiS ) 
and(M*,M* s ) are linearly independent in H~^(dtt,M. 2 ) x H~i(d£l,R 2 ). Let 
Wi, w* , % = 1,2, 6e solutions to the Neumann problem fl4.4ip - fl4.45l) . with 
D = D it corresponding to boundary data M , M* respectively. If 

Wi = w 2 , wi. n = w 2>n , on T, (4.46) 

w{ = w* 2 , wl :Tl = w^ n , on T, (4.47) 

then 

D 1 = D 2 . (4.48) 

Next, let us consider the case of a plate whose boundary is composed 
by an accessible portion T and by an unknown inaccessible portion /, to be 
determined. More precisely, let T, / be two closed, nonempty sub-arcs of the 
boundary dfl such that 

ru/ = ffl, rn/ = {g ;J R}, (4.49) 

where Q, R are two distinct points of d£l. On the assigned couple field M 
let us require the following assumptions: 

MGL 2 (r,R 2 ), (M n ,M r>s )^0, (4.50) 

J % = 0, i = 1,2. (4.51) 

Under the same assumptions made in Subsection 14. II for the plate tensor and 
the domain Q, the transversal displacement w G H 2 (Q) satisfies the following 



Neumann problem 

' = 0, in n, (4.52) 

M a pn a nfs = M n , on T, (4.53) 

< M a p t pna + (M a pnpTa),s=-(^r)„, on T, (4.54) 

M a pn a n p = 0, on /, (4.55) 

^ Map^ria + (M a pn P T a ), s = 0, on I. (4.56) 



Concerning the inverse problem of the determination of the unknown 
boundary portion J, in |M-R3j we have proved the following result. 

Theorem 4.5 (Unique determination of unknown boundaries with one mea- 
surement). Let Qi, Q2 be two simply connected bounded domains in M 2 such 
that dVL i; i = 1,2, are of class C 4 ' 1 . Let dfli — /j U T, i = 1,2, where Ii 
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and T are the inaccessible and the accessible parts of the boundaries dQi, re- 
spectively. Let us assume that f2 x and Q 2 lie on the same side of T and that 
conditions (14.491) are satisfied by both pairs {I±, T} and {I 2 , T}. Let the plate 
tensor P of class C 2,1 (M 2 ) be given by (14 .11) . and satisfying (14.31) (I4.4p and 
the dichotomy condition (I4.5ap or (14.5bl) . Let M G L 2 (r,R 2 ) be a bound- 
ary couple field satisfying conditions (14.501) . (14. 121) . Let Wi G H 2 (Qi) be a 
solution to the Neumann problem (14.521) - (14. 51[) in Q — Qi, i — 1, 2. // 

Wi = w 2 , wi,„ = w 2 , n , on T, (4.57) 

then 

Q 1 = tt 2 . (4.58) 



5 Size estimates for extreme inclusions 

5.1 Formulation of the problem and main results 

Let us assume that the middle plane of the plate Q is a bounded domain in 
R 2 of class C 1 ' 1 with constants po, Mq. In the present section we shall derive 
constructive upper and lower bounds of the area of either a rigid inclusion 
or a cavity in an elastic plate from a single boundary measurement. These 
extreme inclusions will be represented by an open subset D of Q such that 
Q \ D is connected and D is compactly contained in Q; that is, there is a 
number d > such that 

dist(D, dtt) >d p . (5.1) 

In addition, in proving the lower bound for the area of D, we shall introduce 
the following a priori information, which is a way of requiring that D is not 
"too thin". 

Definition 5.1. (Scale Invariant Fatness Condition) Given a domain D hav- 
ing Lipschitz boundary with constants rp and L, where r > 0, we shall say 
that it satisfies the Scale Invariant Fatness Condition with constant Q > if 

diam(D) < Qrp . (5.2) 

Remark 5.2. It is evident that if D satisfies Definition 15. 11 then we have the 
trivial upper and lower estimates 

~ 2 p 2 < area(D) < w 2 Q 2 r 2 p 2 , (5.3) 



'i + VT+L^y 
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where U2 denotes the measure of the unit disk in R 2 . Since we are interested 
in obtaining upper and lower bounds of the area of D when D is unknown, it 
will be necessary to consider also the number r as an unknown parameter, and 
all our estimates will not depend on r. Conversely, the parameters L and Q, 
which are invariant under scaling, will be considered as a priori information 
on the unknown inclusion D. 

For reader's convenience and in order to introduce some useful notation, 
we briefly recall the formulation of the equilibrium problem when the plate 
contains either a rigid inclusion or a cavity, and when the inclusion is absent. 

Let us assume that the plate tensor P G L°°(R 2 , £(M 2 , M 2 )) given by fP]) 
satisfies the symmetry conditions (14. 2p and the strong convexity condition 
([ED . Moreover, let M G H~^(d£l, R 2 ) satisfy fTOTj) . 

If an inclusion D made by rigid material is present, with boundary dD 
of class C 1 ' 1 , then the transversal displacement wr corresponding to the 
assigned couple field M is given as the weak solution wr G H 2 (fl \ D) of the 
boundary value problem 



div(div(PV 2 w H )) = 0, 
(FV 2 w R )n ■ n = -M n , 
div(PV 2 w H ) ■ n + ((¥V 2 w R )n ■ r), 

W R \n G A, 



{M T ), 



dw R 
dn 



dn ' 



in VI \ D, 
on dQ, 
on dQ, 
in D, 

on dD, 



(5.4) 
(5.5) 
(5.6) 
(5.7) 
(5.8) 



coupled with the equilibrium conditions for the rigid inclusion D 



(div(PV\^) • n + (QPV 2 w e R )n -r), s )g- ((fV 2 w R )n ■ n)g, n = 0, 

for every g G A, (5.9) 



where we recall that we have defined w R = w\ n \^ and w l R = w\jj. 

When a cavity is present, then the transversal displacement in Q \ D is 
given as the weak solution wy G H 2 (Q \ D) to the boundary value problem 

div(div(PV 2 w)) = 0, 
(¥V 2 w v )n ■ n = -M n , 

{ div(PV 2 w) • n + ((PV 2 wv> • r), s = (M T ) n 
(FV 2 w v )n ■ n = 0, 

div(PV 2 w v ) • n + {{FV 2 w v )n ■ t), s = 0, 



in Q\D, 


(5.10) 


on dfl, 


(5.11) 


on dQ, 


(5.12) 


on dD, 


(5.13) 


on dD. 


(5.14) 
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Finally, when the inclusion is absent, we shall denote by w G H 2 (£l) the 
corresponding transversal displacement of the plate which will be given as 
the weak solution of the Neumann problem 

div(div(PV 2 u> )) = 0, inQ, / 515 x 

(FV 2 w )n ■ n = -M n , on dQ, (5.16) 

div(PV 2 w ) -n+ ((¥V 2 w )n-r), s = (M T ), S , onffl. (5.17) 

We shall denote by W R , Wy, W the work exerted by the couple field M 
acting on dQ when D is a rigid inclusion, it is a cavity, or it is absent, respec- 
tively. By the weak formulation of the corresponding equilibrium problem it 
turns out that 

W R = - f M T , s w R + M n w Rm = [ _FV 2 w R ■ V 2 w R , (5.18) 
Jan Jn\D 

W v = - I M TjS w v + M n w Vm = [ _PV 2 uy • V 2 w v , (5.19) 
Jdil Jn\D 

W = - f M T , s w + M n w , n = [ PV 2 u> • V 2 w - (5.20) 
Jan Jn 

Note that the works W R , Wy and Wq are well defined since they are invariant 

with respect to the addition of any affine function to the displacement fields 

w R , wy and wo, respectively. In the following, the solutions wy and wo will 

be uniquely determined by imposing the normalization conditions 

[ w v = 0, [ Vuy = 0, (5.21) 

JdD JdD 

I w = 0, I Vw = 0. (5.22) 
Jn Jn 

Concerning the solution w R , we found convenient to normalize it by requiring 
that 

w R = 0, in D. (5.23) 
For a given positive number hi, we denote by D hipo the set 

D hlP0 = {xeD\ dist(x, 3D) > hipo}. (5.24) 

We are now in position to state our size estimates. In the case of a rigid 
inclusion we have the following two theorems. 
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Theorem 5.3. Let Q be a simply connected bounded domain in M. 2 , such that 
dfl is of class C 2,1 with constants po, M , and satisfying (I4.19p . Let D be a 
simply connected open subset ofQ with boundary dD of class C ' , satisfying 
(15.11) . such that Q\D is connected and 

area(D /llPo ) > ^area(D), (5.25) 

for a given positive number h\. Let the plate tensor P given by (14. ip satisfy 
(P , fl44j) . (14^301) and the dichotomy condition. LetMe L 2 (dtt,R 2 ) satisfy 
(14. 26j) - (14.281) . with T satisfying (I4.25p . The following inequality holds 

area(D) < K^ W °~ Wr , (5.26) 

Wo 

where the constant K > only depends on the quantities M , Mi, d , hi, 7, 
Si, M, 5 and F. 

Theorem 5.4. Let Q be a simply connected bounded domain in M. 2 , such that 
dfl is of class C 2,1 with constants po, Mq, and satisfying (I4.19p . Let D be a 
simply connected domain satisfying (15. ip . (I5.2p . such that Q\D is connected 
and the boundary dD is of class C 3,1 with constants rp , L. Let the plate 
tensor P given by (jUJ) satisfy flO}, (H3D and (OD]l . Let M G L 2 (dn,R 2 ) 
satisfy (I4.27p . The following inequality holds 

CPo* ( W °Z Wr ) < area(D), (5.27) 



V W j 
where the function $ is given by 

[0,1) 3 t _.$(t) = (5.28) 

and C > is a constant only depending on Mo, Mi, do, L, Q, 7 and M . 

When D is a cavity, the following bounds hold. 

Theorem 5.5. Let Q be a simply connected bounded domain in M?, such that 
dfl is of class C 2,1 with constants po, M , and satisfying (14. 19 j) . Let D be a 
simply connected open subset ofQ with boundary dD of class C 1,1 , satisfying 
(15.11) . such that Q\D is connected and 

&re&(Dh lPo ) > -area(D), (5.29) 
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for a given positive number hi . Let the plate tensor P given by (14. ip satisfy 
P~Zfl , (l44j) . (ET301 and the dichotomy condition. LetMe L 2 (dtt,R 2 ) satisfy 
(14. 26p - (14.281) . with T satisfying (14.251) . The following inequality holds 

area(D) < Kp 2 V ° , (5.30) 

Wo 

where the constant K > only depends on the quantities M , Mi, d , hi, 7, 
5i, M, 5 and F. 

Theorem 5.6. Let Q be a simply connected bounded domain in M. 2 , such 
that dQ is of class C 1,1 with constants po, Mo, and satisfying (I4.19p . Let 
D be a simply connected domain satisfying (15.11) . (15. 2p . such that Q \ D 
is connected and the boundary dD is of class C 1,1 with constants rp , L. 
Let the plate tensor P given by (14. ip satisfy (14.21) and (I4.4p . and such that 
||P|| C 2,i( R 2) < M' , where M' is a positive parameter. Let M £ L 2 (d£l,M. 2 ) 
satisfy (I4.27p . The following inequality holds 

Cp 2 * ( Wv m W ° ) < areap), (5.31) 
where the function \1/ is given by 

[0,+oo)9t^(t) = T ^, (5.32) 
and C > is a constant only depending on M , Mi, d , L, Q, 7 and M' . 

5.2 Proof of Theorems ET3] and [5T41 

The starting point of the upper and lower estimates of the area of a rigid 
inclusion is the following energy estimate, in which the works Wq and Wr 
are compared. 

Lemma 5.7. Let Q be a simply connected bounded domain in 1R 2 with bound- 
ary dfl of class C 1 ' 1 . Assume that D is a simply connected open set compactly 
contained in Q, with boundary dD of class C 1,1 and such that Q \ D is con- 
nected. Let the plate tensor P £ L°°(fi, £(M 2 , M 2 )) given by (143]) satisfy 
(@2D and ([P]) . Lei M £ i/-^(9fi,M 2 ) safe/y (1427]) . Lei w B £ H 2 (tt\D), 
w £ H 2 (Q) be the solutions to problems (I5.4p - (l5.9p and (15. 15j) - (15.171) . nor- 
malized as above. We have 

[ ¥V 2 w -V 2 wo<Wo-W R = [ M n {w R )wo,n + V{w R )w , (5.33) 

J D JdD 
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where we have denoted by 

M n {w R ) = -{FV 2 w R )n ■ n, (5.34) 

V(w R ) = div(FV 2 w R ) ■ n + ((FV 2 w R )n • r), s (5.35) 

the bending moment and the Kirchhoff shear on dD associated to w R , respec- 
tively. Here, n denotes the exterior unit normal to Q\D. 

Proof. The proof is based on the weak formulation of problems f)5.4l) - fl5.9p 
and fl5.15p - fl5.17p . and can be obtained by adapting the proof of the corre- 
sponding result for a rigid inclusion in an elastic body derived in ( |M-Rlj . 
Lemma 3.1). □ 

Proof of Theorem \5.3[ By (I5.33P and from the strong convexity condition 
ffl~4"j) we have 

/ \V 2 w \ 2 <r\W Q -W R ). (5.36) 

JD 

Estimate f)5.26p can be obtained from the following lower bound for the elastic 
energy associated to w in D 

W 2 w,\ 2 > C^^-Wo, (5.37) 

D Po 

where the constant C > only depends on M , Mi, d , hi, 7, 81, M, 5 
and F. The above estimate was derived in ( |M-R-V e6j. Theorem 3.1) for 
inclusions D satisfying the fatness condition (I5.25P and its proof is based 
on a three sphere inequality for solutions to the plate equation (I5.15P with 
anisotropic elastic coefficients obeying to the dichotomy condition. □ 

In order to prove Theorem 15.41 we need the following Poincare inequalities 
of constructive type. 

For a given positive number r > 0, we denote by D rpo the following set 

D rp0 = {x E R 2 \ < dist(x, D) < rp }. (5.38) 

For D with Lispchitz boundary and u 6 H 1 (D) we define 

Ud = 777T / M ' UdD = T^nT / u - ( 5 ' 39 ) 
\ JJ \ Jd \ au \ JdD 

Proposition 5.8. Let D be a bounded domain in W 1 , n > 2, of Lipschitz 
class with constants rpo, L, satisfying condition (15. 2p with constant Q > 0. 
For every u G H l (D) we have 

[ \u - u D \ 2 < dr 2 p 2 Q [ |Vm| 2 , (5.40) 
Jd Jd 
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\u-u dD y < C 2 rp / |Vw| 2 , (5.41; 
dD Jd 



where C% > 0, C2 > only depend on L and Q. 
Ifue H 1 (D rpo ), then 



[ \u-u dD \ 2 <C 3 rp [ |V«| 2 , (5.42) 

JdD J D r P<) 



where C3 > only depends on L and Q. Moreover, if u G H l (D rp °) and 
u = on dD, then we have 

u 2 < C 4 r 2 p 2 [ \Vu\ 2 , (5.43) 

D r Po J D t po 

where C4 > only depends on L and Q. 

Proof. We refer to [Al-M-Rlj for a proof of the inequalities fl5.40p - fl5.42p and 
for a precise evaluation of the constants Ci, C2, C3 in terms of the scale 
invariant bounds L, Q regarding the regularity and shape of D. 

Inequality (15.4=31) is a consequence of the following result. Let Q be a 
bounded domain in M. n with Lipschitz boundary with constants rp , L, and 
such that diam(f2) < Qrp , Q > 0. Let E be any measurable subset of Q 
with positive Lebesgue measure \E\ > 0. For every v G if x (f2) we have 




v-v E \ 2 < (1 + ^1 Cr 2 p 2 / \Vv\ 2 , (5.44) 

Jn 

where the constant C > only depends on L and Q. The above inequality 
follows from Lemma 2.1 of [A1-M-R4] and from (15.401) applied to the function 
v, see also inequality (3.8) of |A1-M-R4"] . Let us extend the function u G 
H 1 (D rp0 ) into the interior of D by taking u = in D, and let us continue 
to denote by u this extended function, with u G H 1 ^ U D rpo ). Inequality 
flCTjl follows from (f5T44|) by taking v = u, Q = D U D rp0 and E = D. □ 



Proof of Theorem\5.4\ Let g = a + bx + cy be an affine function such that 



the function wq = wq + g satisfies 



w = 0, / Vw = 0. (5.45) 

dD JdD 

The function w G H 2 (Q) is a solution of f l5.15p -( 15.17p . and by (15. 9p . from 
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the right-hand side of (15.331) and by applying Holder's inequality, we have 



W -W R < [ l^\M n (w R )\ 2 j N K B N + 

2 / r \ 2 



\V(w R )\ 2 / \w \ 2 ) =h + h- (5.46) 

dD J \JdD J 

We start by estimating I\. By (I5.4ip and by the definition of wq we have 

/ |wo,n| 2 < Crp / \V 2 w \ 2 < Crp ||V 2 wo||ioo (D) area( J D), (5.47) 

JdD JD 

where the constant C > only depends on L and Q. By the Sobolev em- 
bedding theorem (see, for instance, [Ad] ) . by standard interior regularity 
estimates (see, for example, Theorem 8.3 in [M-R- VeT] ) . by Proposition 15.8} 
by (S3D and by (15T2"0~|) we have 

||VV)|U~(Z)) < ^IKII^(Q) < - ( [ \V 2 w \ 2 ) 2 < -Wj, (5.48) 
Po Po \Jn J Po 

where the constant C > only depends on M , Mi, do, 7 and M. 

The integral J dD \M n (w R )\ 2 can be estimated by using a trace inequality 
(see, for instance, [L-Mj ) and a if 3 -regularity estimate up to the boundary 
dD for w R (see, for example, Lemma 4.2 and Theorem 5.2 of [M-R- Ve5] ) : 

3 2 

JdD r Po i= o JD^T r Po i=0 Jd^po 

(5.49) 

where the constant C > only depends on M , Mi, L, Q, 7 and M. The 
function belongs to H 2 (D rp0 ) and, by (15.231) . w R = and Vw R = on 
dD. Therefore, by applying twice (I5.43p . by (I4.4p and by (I5.18P we have 

2 

y>Pof / \^w R \ 2 < Cr 4 p 4 [ \V 2 w R \ 2 < Cr 4 p 4 W R , (5.50) 

i=0 JD r P0 Jd t po 

where the constant C > only depends on L, Q and 7. By (I5.47p . f)5.48p . 
(J53SD, (E30D we have 

/1 < — W^W^(aiea(D))^ (5.51) 
Po 

where the constant C > only depends on M , Mi, do, L, Q, 7 and M. 
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The control of the term I 2 can be obtained similarly. By a standard 
Poincare inequality, by (I5.4ip and by ( I5.48P we have 

/ \w \ 2 < CV 2 po / \wo, s \ 2 < Cr 3 p 3 [ \V 2 w \ 2 < 

JdD JdD JD 

< Cr 3 pol|V 2 Wo||ioo (D) area(D) < Cr 3 'p W area(.D), (5.52) 

where the constant C > only depends on M , Mi, d , L, Q, 7 and M. 
Concerning the integral J dD \V(wr)\ 2 , by using a trace inequality, a ir- 
regularity estimate up to the boundary dD for wr (see, for example, Lemma 
4.3 and Theorem 5.3 in |M-R-Ve5] ) and by (I5.50p . we have 

/ i^)i 2 <^X>of / iwi 2 < 

JdD r Po i=Q JD^T 

^ "Fa / l V V ^ -TaW R , (5.53) 

where the constant C > only depends on M , M 1; L, Q, 7 and M. Then, 
by (I5.52p and ( I5.53P we have 

h < —W^Wl (area(D))*, (5.54) 
Po 

where the constant C > only depends on M , Mi, do, L, Q, 7 and M. By 
fl53Tj) and (15341) . the inequality ( l5T2Tj) follows. □ 

5.3 Proof of Theorems 15.51 and 15.61 

As in the proof of the upper and lower estimates of the area of a rigid inclu- 
sion, we need to compare the works Wo and Wy- The analogue of Lemma 
15.71 is the following result. 

Lemma 5.9. Let Q be a simply connected bounded domain in IR 2 with bound- 
ary dfl of class C 1,1 . Assume that D is a simply connected open set compactly 
contained in Q, with boundary dD of class C 1,1 and such that Q \ D is con- 
nected. Let the plate tensor P E L°°(Q, £(M 2 , M 2 )) given by (14.1 j) satisfy 
(TC2D and f|Q|l . Let M e (dVL, R 2 ) sate/?/ (TC2ZD. Let w v e H 2 (tt\D), 
w e H 2 (Q) be the solutions to problems (I5TTU]) - fl5TT4"P and floTTBD - flSTTTj) . 
normalized as above. We have 

/ PVV, • V 2 w < Wy - Wo = / M n (w )w v ,n + V{w )w v , (5.55) 

J D JdD 

where the functions M n (wo), V{wq) are defined as in (I5.34p and (I5.35p . and 
n denotes the exterior unit normal to Q\D. 
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Proof. As for Lemma 15.7} the proof is based on the weak formulation of 
problems (I5.10p - fl5.14p and fl5.15p - fl5.17p . and can be obtained following the 
same guidelines of the corresponding result for a cavity in an elastic body 
derived in QM-R1] . Lemma 3.5). However, here we simplify the approach 
presented in |M-R1] without extending the function Wy in the interior of 
D. □ 

Proof of Theorem \5.5[ The proof follows from the left-hand side of (I5.55P by 
using the same arguments as in the proof of Theorem 15.31 □ 

Proof of Theorem \5.b\ From the right-hand side of (15.551) and by applying 
Holder's inequality, we have 



W v -W <^J^\M n (w )\ 2 j y] dD \ w vAj + 

+ (JdD |y(W ° )|2 ) 2 (J BD K ' 2 ) 2 =Jl + J * (5 - 56) 

Let us estimate the integral J^- By the Sobolev embedding theorem (see, 
for instance, |Ad] ) , by standard interior regularity estimates for w (see, for 
example, Lemma 1 in |M-R-Ve2j ). by (I5.40p . by f O|) and by (I5.20p we have 

/ \V(w )\ 2 < C\\V 3 w \\ 2 LxiD) \dD\ < ^\\w \\ 2 H2{n) \dD\ < 

JdD Po 

<-A \V 2 w Q \ 2 \dD\ < ^W \dD\, (5.57) 
Po Jn Po 

where the constant C > only depends on M , Mi, d , 7 and M' . By 
( |A1-Rj . Lemma 2.8) we have 

\dD\ < C aTe < D \ (5.58) 
^Po 

where the constant C > only depends on L. 

To control the integral j dD \ wv\ 2 we use a standard Poincare inequality on 
dD, Proposition 15.81 inequality (I5.42p . the strong convexity condition (14. 4p 
for P and the definition of Wy, that is 

/ \w v \ 2 <Cr 2 p 2 [ \w v J 2 < Cr 3 p 3 [ \V 2 w v \ 2 < 

JdD JdD Jd t po 

< Cr 3 pl [ _FV 2 w v ■ V 2 uy = Cr 3 p 3 W v , (5.59) 
Jn\D 
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where the constant C > only depends on L, Q and 7. Then, by (|5.57p . 
(I5.58P and (15.59 p . and since, trivially, r < K, where K > is a constant only 
depending on L and Mi, we have 

■h < — W^Wv(area(D))*, (5.60) 
Po 

where the constant C > only depends on M , M 1; do, L, Q, 7 and M 1 . By 
using similar arguments we can also find the analogous bound for Ji, and 
the thesis follows. □ 
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